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Bright features ('knots') in extragalactic jets are commonly associated with re- 
heating of charged particles in regions with shocks or strong compression. We 
find a natural appearance of nozzles in simulations of toroidally magnetized jets, 
and propose this as a generic mechanism of knotted structures such as in the jet 
of 0800+608. Nozzles provide correlated enhancement of emission in radio and 

^ \ optical, such as discovered in HST observations of 3C273. The simulations re- 
veal a rich structure of knots, including a nose cone with a propagating knotted 
structure and a hot-spot as the most extended brightness feature. Unexpected 
is the observed separation of the hot-spot from a composite Mach disk closer to 
O ' the source. 

Enhanced brightness features in the form of localized 'knots' are perhaps the most signifi- 

.^J ' cant morphological aspects of astrophysical jets. These features are produced with relativis- 
tic motion in both extragalactic sources (e.g. 3C273 1 and M87 2 ) and recently discovered 

d ; 

galactic sources (GRO J1655-40 3 ' 4 and GRS 1915-105 4 ' 5 ). On the larger scale, the most 
extended brightness feature (the 'hot-spot') has been explained as enhanced synchrotron 
radiation behind a terminal shock 6 ' 7 . Increased resolution and sensitivity of radio observa- 
tions shows bright features over the entire length of some jets, from close to the source up to 
the terminal hot-spot, which can be interpreted likewise 8,9 ' 10 ' 11 . Knotted structures appear 
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with extraordinary dynamic range in brightness and patterns, including knots as the first 
detectable features at finite angular distance away from the source (leaving an apparent gap 
between source and jet, as e.g. in M84 12 ), and regular or quasi-periodic knotted structures 
which characterize the large scale morphology (e.g. as in 3C 111 13 , QSO 0800+608 14 ' 15 (see 
Figure 1), Cyg A 16 ). If knot motion indeed relates to bulk motion, BL Lacertae objects, 
radio quasars and radio galaxies emit plasma with appreciable Lorentz factors 17 . Relativistic 
effects introduce beaming (cf. Ref 11), which allows for the unification scheme for compact 
radio loud sources. 

In this Letter, the nature of knots is considered in the large scale morphology of jets 
modeled by the theory of time-dependent, fully relativistic ideal magneto- hydrodynamics 
(MHD) of a polytropic, perfect fluid as the result of boundary conditions associated with 
the 'standard model' of active galactic nuclei (a rotating black hole and accretion disk 24 ' 25 ; 
see also Refs 26 and 21 for a discussion). The simulations use the equations of MHD in 
divergence form 18 " 23 . 

Quasi-periodic structures have a natural appearance in simulations of toroidally magne- 
tized relativistic jets, which are shed off by an oscillation in the hot-spot in the head of the 
jet subsequent to a stagnation point near the source (Figure 2 21 ). Figure 2 shows a stag- 
nation point (first, left pressure maximum), and three subsequent local pressure maxima 
with Mach number M — 1, which are sonic nozzles with pressure and density contrasts 
4.39 ±L37 and 2.67 ±a57 , respectively. The nozzles constitute funnels of nearly unmagnetized 
flow, pinched by the magnetized flow at larger radius; since the toroidal magnetic field van- 
ishes on axis, it is dynamically appreciable only some distance away from the axis and the 
flow near-axis is to leading order purely hydrodynamic. The basic scaling of the nozzles can 
be derived by considering a tube of streamlines through a characteristic cross section, A s , at 
the sonic point, and its (for supersonic flow wider) cross section, A, in between two nozzles. 
The density contrast y as a function of the ratio of widths a = is parametrized by 
the sound speed, a s , at the sonic point. The nonrelativistic (a s — > 0) and ultra-relativistic 
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(a s — > \/7 ~ 1) limits thus provide upper and lower bounds, respectively, for f- close to the 
axis, which for large a are asymptotically given by 



Figure 3 shows the result of a supersonic (M = 1.67), low density relativistic (T = 2.46) 
jet with out of radial force-balance boundary conditions due to a toroidal magnetic field as 
the "spin-off" in the energy extraction process by poloidal field lines in the source. The jet 
propagates into an unmagnetized outer medium at rest and in hydrostatic pressure balance 
with the jet. The central jet flow partially thermalizes in a Mach disk, which consists of a 
shock of small radius, which is by-passed by jet flow which shocks in an annular shock further 
downstream. In between these two shocks resides a much weaker and larger shock wave, 
the combination of which with the small Mach disk we shall refer to as the composite Mach 
disk. The shocked by-pass flow subsequently bifurcates into a cocoon and a recollimated 
nose cone. The cocoon exerts a pinch on the boundary of the jet in the neighborhood of the 
small Mach disk, which initiates the formation of a weakly shocked nozzle; a feedback loop 
is thus obtained (see Ref 27 for a feedback loop in hydrodynamical jets). Further refinement 
of the computations are needed to detail this process. At the 'root' of the nose cone resides 
a stagnation point (as in the nonrelativistic simulations of Refs 28 and 29), which sets 
the stage for the formation of nozzles as in Figure 2 in the nose cone ahead of the Mach 
disk. A hot-spot results as the most extended enhanced brightness feature, separated off 
the composite Mach disk and flow bifurcation structure. Exploration in parameter space is 
needed to obtain critical values at which this morphology is typical. Because the nozzles form 
on axis by the internal dynamics of the jet, the formation process is expected to be rather 
insensitive to (the interaction with) the environment, and hence the process is expected to 
persist in three-dimensional simulations (i.e. with no imposition of cylindrical symmetry). 
This simulation illustrates that we can learn from the extended morphology (up to the 
terminal hot-spot) about the characteristics of the flow at the source (the Mach number and 
radial distribution of toroidal magnetic field), which may serve as a new approach towards 
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probing the energy extraction process in the core of active galactic nuclei. 

If the knots of 0800+608 turn out to be (i) standing still, our results predict the jet to 
be transonic (with a stagnation point close to the source; Figure 2); (ii) moving, our results 
predict the jet to be supersonic (with the stagnation point moving ahead of the Mach-disk; 
Figure 3). The asymmetry in the morphology of 0800+608 has been attributed by Jackson 
et a/. 15 to a wind interacting disruptively with the counter jet, which favors (i) if the jets are 
symmetric at the source. In alternative (ii), the recollimation associated with the stagnation 
point may denote a transition to constant width in an environment of decreasing pressure 
(knot 3 in Figure 1). The weak extended structure beyond knot 8 then corresponds to part 
of the thin cocoon around the nose cone in Figure 3. Jackson et a/. 15 infers the presence of 
a cocoon between knot 3 and knot 8 from edge-brightening due to jet-cocoon interaction, 
which would need further refinement in our model to be brought about in emissivity. Our 
results support the suggestion of Jackson et a/. 15 of the knots in 0800+608 being associated 
with an intrinsic, stable and periodic magnetic pinch, as opposed to being due to an unstable 
interaction with the environment. 

The adiabatic process associated with the sonic knots finds a place in the observed 
correlations between radio and (HST-)optical synchrotron emission in the recent observations 
on 3C273 30 ' 31 . The relativistic simulations reported here also serve as a step towards fully 
three-dimensional simulations with helical magnetic fields, as required in the unification 
scheme in compact radio loud sources. 
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FIGURE 1. A 6cm (0.35-arcsec circular beam) total intensity VLA map of 0800+608 
(,2=0.689, L=140kpc with H=50kms~ 1 Mpc~ 1 ) by Jackson et. al. (1990). Contours are from 
0.05-6.4 mJy beam -1 (Courtesy of Blackwell Science Ltd). 

FIGURE 2. Coordinate distributions of pressure (A) and magnetic field strength (B) of a 
fully developed axisymmetric relativistic jet (v < 0.73c, where c is the velocity of light), pro- 
duced by a Newtonian sphere with boundary conditions (H,P,r) = (0.72 sin 9, 0.100, 0.200) 
(in spherical coordinates with 9 = on the axis of symmetry). The fluid has 7 = §, and 
is released from the surface of the sphere by diffusion. The Mach number M = - on-axis 
(C) (sound speed a = {-pj}^) shows the formation of sonic nozzles; solid and dashed lines 
correspond to two consecutive times, 15% apart, with local Mach numbers at the three 
sonic points (in between brackets give values at the later time). The jet shows a stagnation 
point (first, left knot), and subsequently three sonic (M = 1) nozzles in various degrees 
of maturation, shed off by an oscillation in the hot-spot in the head of the jet. Schematic 
diagram (D) of a tube of streamlines in a sonic (M = 1) nozzle with cross section A s at 
the sonic point and A at maximal expansion. The flow is described by Bernoulli's equa- 
tion, H = f*T = f*T s , and continuity, $ = ru z A = r r u z A s , where /* = / + kr 2 and the 
subscript s refers to the sonic section. Because k = £ = O(a), the near- axis flow is to 
leading order a hydrodynamical flow, and /* ~ / facilitates the relationship of the density 
contrast ^ as a function of a = {-^-)^ ■ Using u b ~ (cosh A, sinh A, 0, 0) near axis, the non- 
relativistic limit a s -> gives a (M) = (^ T )^ T )M~2(1 + ~ (^) J ~M~ 

as M — > oo, where M = (^-) 5_ ^2, which serves as an upper bound on y. The rela- 
tivistic case retains the sound speed, tanhA s = a s , at the sonic point as a parameter in 
oti(^;a s ) = (^)~5{-j|cosh 2 \ s - l}-i s inh5 A s ~ r'^r 1 ^ - l)i as a s -> (7 - 1)5, where 

s 2 2 

7 = 1 + ^zf (^j) 7 " 1 with f s — (1 — t^tj-)" 1 . The computations are in cylindrical coordinates 
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(t, a, (ft, z), using the covariant and constraint-free divergence formulation of relativistic MHD 
(van Putten 1991-1995), and are performed on the Intel Paragon parallel computer at SDSC, 
San Diego. 

FIGURE 3. Coordinate distributions of velocity magnitude, pressure, rest mass density 
and magnetic field strength of fully developed axisymmetric, toroidally magnetized jet at 
time t/aj e t = 34.47, where aj et is the radius of the jet (in arbitrary units). The colors 
vary linearly from blue to red. The jet aperture has boundary conditions (r, M, H, P, r) = 
(2.46, 1.67, 0.46<7 cos -f— t \, 0.100, 0.200) which are out of radial force-balance, and the en- 
vironment is unmagnetized and initially at rest with (P, r) = (0.100,1.00). The aperture 
boundary conditions at a — <jj et are smoothed through multiplication by a tangent hyper- 
bolic (this incidentally gives a thin layer of reversed toroidal magnetic field). The on-axis 
distributions of pressure and rest-mass density of the jet vary by a factor of 190 and 45, 
respectively. The solution is characterized by a hot-spot of postshock material at the head 
of the nose cone (as in Refs 6 and 7), which has separated off from the composite Mach 
disk. The oscillation in the terminal shock produces a propagating (v = 0.162 ±0 04 c), super- 
sonic nozzle (M = 1.28 ±003 in the co- moving frame) with pressure contrast 6.14 and rest 
mass density contrast 3.38. A repeat of the oscillation in the terminal shock is observed 
at t/(jj e t = 39.06, which produces a second nozzle (not shown). The computations use the 
same method of that for Figure 2, and are performed on the IBM SP2 parallel computer at 
the Cornell Theory Center. 
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